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ABSTRACT 

Gravitational lensing predicts an enhancement of the density of bright, distant QSOs 
around foreground galaxies. We measure this QSO-galaxy correlation Wqg for two 
complete samples of radio-loud quasars, the southern IJy and Half-Jansky samples. 
The existence of a positive correlation between z ~ 1 quasars and z « 0.15 galaxies 
is confirmed at a p = 99.0% significance level (> 99.9% if previous measurements on 
the northern hemisphere are included). A comparison with the results obtained for 
incomplete quasar catalogs (e.g. the Veron-Cetty and Veron compilation) suggests the 
existence of an 'identification bias', which spuriously increases the estimated ampli- 
tude of Wqg for incomplete samples. This effect may explain many of the surprisingly 
strong quasar-galaxy associations found in the literature. Nevertheless, the value of 
Wqg that we measure in our complete catalogs is still considerably higher than the 
predictions from weak lensing. We consider two effects which could help to explain 
this discrepancy: galactic dust extinction and strong lensing. 
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1 INTRODUCTION 

Canizares (1981) showed that if gravitational lensing effects 
are considered, an enhancement in the density of QSOs close 
to the position of foreground galaxies is expected (see also 
Keel 1982; Peacock 1982). This effect, known as 'magnifica- 
tion bias' can be characterized by the overdensity q — jj."^^ 
(Narayan 1989), where jj, is the magnification and a is the 
logarithmic slope of the number counts distribution. De- 
pending on the slope q, an excess (g > 1) or even a de- 
fect {q < 1) of background sources around the lenses will 
be observed. Galaxies and clusters trace the matter fiuctua- 
tions (up to a bias factor b), and therefore, we may observe 
positive (q > 1), null (a = 1) or negative (a < 1) statis- 
tical associations of these foreground objects with distant, 
background sources. 

The evidence for association between high-redshift 
AGNs and foreground galaxies has steadily accumulated 
during the years (see Schneider, Ehlers & Falco 1992 for 
a detailed review). However, the results are sometimes ap- 
parently contradictory and quite often difficult to accom- 
modate within the gravitational lensing framework. For in- 
stance, the strength of associations depends on the AGN 
type and in particular, the studies performed with radio- 
loud AGN samples, or with heterogeneous samples, usually 
extracted from early versions of the Veron-Cetty & Veron 
or the Hewitt & Burbidge catalogs, (which contain a high 
proportion of radio-loud QSOs) almost routinely find signif- 
icant excesses of foreground objects around the quasar posi- 



tions (Tyson 1986; Fugmann 1988,1990; Hammer & Le Fevre 
1990; Hintzen et al. 1991; Drinkwater et al. 1992; Thomas 
et al. 1995; Bartelmann & Schneider 1993b, 1994; Bartsch, 
Schneider, & Bartelmann 1997; Seitz & Schneider; Bemtez et 
al. 1995; Bem'tez & Martinez-Gonzalez 1995, 1997(BMG97); 
Bem'tez, Martinez-Gonzalez & Martm-Mirones 1997; Nor- 
man & Williams 1999; Norman & Impey 1999). Although 
these results are qualitatively in agreement with the magni- 
fication bias effect, in most cases the amplitude of the corre- 
lation is much higher than that expected from gravitational 
lensing models. 

On the other hand, the studies carried out with opti- 
cally selected catalogs are less consistent. Positive correla- 
tions have been found by Webster et al. (1988), Rodrigues- 
Williams and Hogan(1991) and Williams & Irwin (1998). 
In all these cases, the amplitude of the correlations is more 
than an order of magnitude stronger than the predictions 
from the magnification bias effect. On the other hand, Boyle 
et al. (1988); Romani & Maoz (1992); BMG97; Ferreras, 
Bem'tez & Martinez-Gonzalez (1997) and Groom & Shanks 
(1999) found null or even negative correlations, which in 
some cases are expected from the magnification bias efi'ect 
due to the shallowness of the QSO number counts. However, 
there are certain instances in which the differences in the 
QSO-galaxy correlation Wqg between radio-loud and radio- 
quiet QSOs cannot be explained by the lensing hypothesis: 
BMG97 found a negative correlation for the optically se- 
lected Large Bright Quasar Survey (LBQS) catalog (Hewett, 
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Foltz & Chafec 1995), which has a steep slope a ^ 2. The 
most hkely explanation for this is that the LBQS catalog is 
affected by a selection bias which hinders the detection of 
high-z QSOs in regions of high projected galaxy density (see 
Roinani & Maoz 1992, Maoz 1995). Further evidence was 
obtained by Ferreras, Beni'tez & Martmez-Gonzalez (1997) 
which found that CFHT/MMT QSOs (Crampton, Cowley 
& Hartwick 1989) were strongly amticorrelated with low red- 
shift galaxies, contrary to what was expected from the mag- 
nification bias effect. The amplitude of the anticorrclation is 
stronger for lower redshift QSOs which practically excludes 
dust in foreground galaxies as its cause and strongly suggests 
a dependence on the strength of the emission lines used to 
identify the QSO. From these results it may be concluded 
that until the selection effects operating on optically selected 
QSO samples axe not better understood, one should be wary 
of the gravitational lensing inferences obtained from them. 

As Bartelmann & Schneider (1993a) pointed out, the 
scale of some of the positive correlations reported above 
(several arcmin) is difficult to explain considering lensing 
by isolated galaxies or microlensing, and apparently has to 
be caused by the large-scale structure. Bartelmann (1995) 
showed that in the weak lensing regime 



Wgg{e) = b{a - l)C^s 



(1) 



where C^s is the correlation between magnification and mat- 
ter density contrast C^jf SfiSm >, & is the biasing factor 
b = 5g/Sm, where Sg and 5m are respectively the galaxy and 
dark matter overdensities. Sanz et al. (1997) and Dolag & 
Bartelmann (1997) have calculated C^a for several cosmo- 
logical models taking into account the non linear evolution 
of the power spectrum of density perturbations. In Bemtez 
& Sanz (1999) it was shown that the expected value of Wqg 
can also be estimated as 



Q{Zb, Zf)Qmb ^Wg 



(2) 



where Zf and Zb are the typical redshifts of respectively the 
galaxy and quasar samples. For low redshift galaxies and rel- 
atively high-z QSOs like the ones considered in this paper, 
the factor Q is approximately independent of the cosmolog- 
ical model. 

Any radio-loud quasar is detected in the optical, with 
e.g. magnitude B, and in radio with flux S. For a sample 
with independent optical and radio fluxes, the effective slope 
Qe (Borgeest et al. 1991) takes the form 



ae{S,B) 



dlnN{>S) 2.5dlogAf(< B) 



d\nS 



dB 



(3) 



where N{> S) and A'"(< B) are the cumulative number 
counts distributions. This is known as the double magnifi- 
cation bias effect. It holds if the optical and radio fluxes are 
independent, and also if the source sizes are considerable 
smaller than the lens in both bands, what ensures that both 
fluxes are equally magnified. Both conditions arc reason- 
ably fulfilled for QSO-galaxy associations: radio and optical 
fluxes arc practically uncorrelated for radio quasars (see e.g. 
Fig. 10 of Drinkwater et al. 1997) and the large scale struc- 
ture weak lensing field is smooth enough to magnify similarly 
the sizes of the optical and radio emitting regions. 

Therefore, to compare the model predictions with the 
observations, it is necessary to know in detail the unper- 
turbed number counts distribution of the quasar sample. 



However, practically all the radio-loud samples mentioned 
above are incomplete (e.g. BMG97) or have heterogeneous 
photometrical information in the optical (e.g. Bemtez & 
Martmez-Gonzalez 1995). Here we intend to remedy this 
situation by measuring Wqg using two complete radioloud 
quasar samples. 

The outline of the paper is the following. In §2 we de- 
scribe the galaxy and QSO samples used to estimate Wqg. 
§3 deals with the statistical analysis. §4 presents two theo- 
retical estimations of Wqg and compares them with the ob- 
servational results. §5 discusses the puzzling results of the 
previous section and §6 summarizes our main results and 
conclusions. 



2 THE DATA 

Our galaxy sample is obtained from the ROE/NRL COS- 
MOS/UKST Southern Sky catalog (see Yentis et al. 1992 
and references therein), which contains the objects de- 
tected in the COSMOS scans of the glass copies of the 
ESO/SERC blue survey plates. The UKST survey is orga- 
nized in 6 X 6 deg^ fields on a 5 deg grid and covers the zone 
with declination —90 deg < S < Odeg and galactic latitude 
|6| > 10 deg. The catalog supplies several parameters for 
each detected object, including the centroid in both sky and 
plate coordinates, Bj magnitude and the star-galaxy image 
classification down to a limiting magnitude of Bj w 21. 
Drinkwater et al. (1997) quote a calibration accuracy of 
about ±0.5mag for the COSMOS Bj magnitudes. As in 
BMG97, we include in our sample galaxies brighter than 
Bj < 20.5 and within 0.2' < 6 < 15' of the quasars selected 
as we describe below. The median redshift of the Bj galaxies 
is Zg « 0.15. 

In BMG97 we studied the quasar-galaxy correlation for 
a sample of 144 z > 0.3 PKS quasars with a 2.7GHz fiux 
82.7 > 0.5Jy, extracted from Veron-Cetty & Veron (1996). 
This sample was shown to be strongly correlated with fore- 
ground galaxies, but its incompleteness precluded a detailed 
comparison with theory. Thus, we now consider two well de- 
fined and practically complete radio loud quasar catalogs 
to assess the incompleteness of the BMG97 sample and its 
effect on the estimation of Wqgi the IJy catalog (Stickel, 
Meisenheimer & Kiihr 1994; Stickel et al. 1996) which in- 
cludes radiosources with a 5GHz flux S5 > IJy, and the 
Parkes Half-Jansky Flat-Spectrum sample (Drinkwater et 
al. 1997) which contains flat-spectrum sources detected at 
2.7GHz and 5GHz, with S2.7 > 0.5Jy, spectral index be- 
tween 2.7GHz and 5GHz Q2. 7/5.0 > ~0.5, galactic latitude 
|fe| > 20 deg and declination —45 deg < S < 10 deg. 

Both the IJy and Half-Jansky catalogs are almost 100% 
identified and have spectroscopical redshifts for > 90% of 
the sources. To try to work with samples as similar as pos- 
sible, we have extracted from the three catalogs those ob- 
jects classified as QSOs. VC classify as QSOs those objects 
with absolute magnitudes brighter than B = —23 and the 
IJy catalog uses the same classifications as the original dis- 
covery papers. For the HJ we have selected those objects 
classified as "stellar" , excluding the BL LACs to get a more 
homogeneous sample. Drinkwater & Schmidt 1996 showed 
that the original Parkes catalog classification was made in 
a non-uniform way, resulting in false large-scale correlations 
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among quasars drawn from that catalog. The HJ classifi- 
cation seems to be free of this problem (Drinkwater ct al. 
1997) and even if the IJy presents a similar defect, it would 
not affect the measurement of Wqg, since we normalize the 
galaxy density locally in a 15' circle around each quasar. It 
should also be kept in mind that the optical classification of 
radiosources may be affected by crowding effects (an exam- 
ple are the 'merged' objects in Drinkwater et al. 1997), and 
that this could tend to eliminate from the sample QSOs in 
fields with higher than average galaxy density, biasing low 
the measured value of Wqg . 

Our analysis will be performed within basically the 
same area defined by Drinkwater ct al. (1997), except for 
three additional constraints imposed by the characteristics 
of the COSMOS/UKST galaxy catalog: 

i) Declinations — 45deg < 5 < 3deg, as the COS- 
MOS/UKST plates only reach up to <5 < 3 deg. 

ii) Galactic latitude |&| > 30 deg. It was shown in 
BMG97 that the "galaxy" density in COSMOS fields is sig- 
nificantly anticorrelated with galaxy latitude, clearly due to 
an increase in the numbers of stars misclassified as "galax- 
ies" when we approach the galactic disk. In fact, COSMOS 
fields with |6| < 30 deg have on average 75% more objects 
classified as "galaxies" than the |b| > 30 deg fields. Although 
excluding these fields does not afi'ect strongly the results 
(BMG97), they would dilute the galaxy excess and bias low 
the amplitude of Wqg . 

ifi) \Ax\,\Ay\ < 2.25 deg, ^Ax^ + Aj/2 < 2.5 deg. As in 
BMG97, these cutoffs avoid the outer regions of the plates, 
with worse image and photometric qualities. We also exclude 
a few fields from the sample because they are affected by 
meteorite traces. 

There are also other three general constraints based on 
the quasar characteristics: 

iv) A faint threshold on the quasar Bj magnitude, 
Bj < 20.5, the same as for our galaxies. We also sot an up- 
per threshold of Bj > 15, brighter of which the photographic 
magnitudes are not reliable. We use only COSMOS Bj mag- 
nitudes obtained from Drinkwater et al. (1997), or directly 
from the ROE/NRL COSMOS/UKST Southern Sky cata- 
logue. This ensures photometrical bandpass homogeneity. At 
the chosen magnitude limit, it can be resisonably expected 
that any radio source candidate is bright enough to be de- 
tected, correctly identified on a photographic plate and have 
its redshift determined spectroscopically. Up to this limit, 
the quasar catalogs should be practically complete. 

v) A lower redshift cutoff, z > 0.3. Only 5% of the COS- 
MOS/UKST galaxies with Bj < 20.5 have z > 0.3, so we 
exclude the possibility of "intrinsic" quasar-galaxy corre- 
lations contaminating the result. Setting a higher redshift 
cutoff does not affect significantly the results. 

vi) A 2.7GHz flux cutoff: S2.7 > 0.5 Jy. This is the flux 
limit of the Half-Jansky catalog. The 52.7 fluxes are obtained 
from Drinkwater et al. (1997) and from the Veron-Cetty & 
Veron compilation. Objects from the IJy catalog also have 
S5 > 1 ,Jy. We do not set any constrains on the steepness 
of the radio spectra for the Veron-Cetty & Veron and IJy 
samples. The Half-Jansky sample has 02.7/5.0 > —0.5. 

The Veron-Cetty & Veron (VC) sample is very similar 
to the PKS sample used in BMG97 except for the exclusion 
of 23 quasars with |6| < 30 deg and Bj < 20.5, and the 
inclusion of 12 objects which belong to the Veron-Cetty & 




Figure 1. Redshift histograms for the VC (solid line), HJ (daslied 
line) and IJy (dotted Une) samples. The vertical axis show the 
fraction of quasars in each redshift bin. 



Veron (1997) catalog and comply with criteria i-vi) but were 
not included in BMG97 because their names do not start 
with "PKS", or have S2.7 exactly equal to 0.5Jy (the flux 
limit in BMG97 was 52.7 > 0.5Jy and here it is S2.7 > 
0.5 Jy). The VC sample thus defined contains all the quasars 
in the IJy sample and 91.5% of the Half-Jansky (HJ) sample 
(86 out of 94). 

There are 35 VC quasars which do not belong to the HJ 
or IJy sample. They are steep-spectrum quasars fainter than 
the IJy catalog > 1 Jy limit, and are slightly brighter in 
the optical (< Bj >= 17.50) than the HJ sample. On the 
other hand, there are 8 HJ quasars which are not included 
in the VC sample. They are considerable fainter than the 
rest of the HJ sample, with Bj > 18.6 and 6*2.7 < 0.8. 

The redshift distributions of the HJ, VC and IJy sam- 
ples are shown in Fig 1. The distributions look rather simi- 
lar, although the HJ quasars apparently tend to have slightly 
higher redshifts than the IJy ones. However, a Kolmogorov- 
Smirnov test shows that the difference between both distri- 
butions has only a 69% significance level. As Drinkwater et 
al. (1997) showed, the redshift distribution of the HJ cata- 
log is very similar to that of optically selected samples, as 
the LBQS (Hewett, Foltz & Chafee 1995). Fig 2 displays 
the Bj histograms of the HJ, VC and IJy samples. Al- 
though VC and IJy quasars are slightly brighter than HJ 
ones, the difference is again not statistically significant. The 
cumulative number distributions N{< Bj) of the z > 0.3 
quasars in the HJ and IJy samples are shown in Fig. 3. 
They are much shallower than the corresponding distribu- 
tion for optically-selected QSOs (Hawkins & Veron 1995). 
This fact was also noticed in BMG97 but since we were us- 
ing a subsample of an incomplete catalog, it was not possible 
to distinguish between the effect of incompleteness and the 
intrinsic differences in the luminosity functions. The number 
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Figure 2. Bj histograms for the VC, HJ and IJy samples. The 
vertical axis show the fraction of quasars in each bin. 

counts A'^(< Bj) for the quasars are well fitted by a law of 

the form 

logiV(< Bj) = log Nb + aiAB J - azAB J, ABj < 

logAf(< Bj) = logA^s +aiAB.7, ABj > 

where ABj = Bj — 18.75 and Nb is an arbitrary normal- 
ization. For the IJy sample ai = 0.19, 02 = 0.66 , whereas 
for the HJ sample, ai = 0.11, 02 = 0.77. The above fits are 
plotted in Fig. 3. 

The cumulative flux distribution of the z > 0.3 quasars 
in the HJ and IJy catalogs axe shown in Figure 4. For the 
HJ sample 

log N{> 82.7) = log TVo - 1.46 log 52.7 - 0.36(log 82.7)^ (4) 
and for the IJy sample 

log N{> S5) = log No - 1.82 log Ss (5) 

where A'o is an arbitrary normalization. 

Therefore, at the flux limits of the HJ sample, Bj = 
20.5 and 52.7 = 0.5Jy, we have that (oe - 1)hj = 0.18 + 
1.24 - 1 = 0.42. For the IJy sample, (q^ - l)ijy = 0.11 + 
1.82 - 1 = 0.93. The real value of Oe - 1 is therefore 2-4 
times smaller than the values so far considered in the esti- 
mation of Wqg and q for the PKS sample (BMG97, Dolag & 
Bartelmann 1997), Qe - 1 = 1.8 - 2.2. 



3 STATISTICAL ANALYSIS 

Spearman's rank order test has often been used to study 
the statistical significance of QSO-galaxy associations. For 
instance, in the implementation of Bartelmann & Schnei- 
der (1993a, 1994), all the individual galaxy fields around 
the quasars are merged into a single 'superfield', which is 



Figure 3. Cumulative number counts-magnitude distribution for 
the HJ and IJy samples. The solid lines are the least squares fits 
described in Sec 2 



Table 1. Results of the weighted average test 



Sample 


Nq 






Pw 


VC 


133 


12662 


1.01375 


99.55% 


IJy 


54 


5339 


1.01537 


97.09% 


HJ 


94 


9243 


1.01178 


97.33% 


HJ + IJy 


106 


10415 


1.01200 


97.84% 


VC-l-HJ-l-lJy 


141 


13577 


1.01134 


98.57% 


VC "faint" 


8 


553 


1.03805 


93.43% 


HJ "faint" 


14 


1309 


0.99953 


50.37% 


"non-VC" HJ 


8 


915 


0.97132 


6.72% 



Nq is the number of QSOs in each sample, Ng the number of 
galaxies for which the Wqg is measured, r-w is the value of the 
normalized correlation coefficient (sec text) and pu, the corre- 
sponding significance associated to the value of Vw ■ 

subsequently divided into Nuns annular intervals of equal 
surface. The rank order test determines whether the num- 
ber of galaxies in each bin rii, (i = 1, Ntins) is anticorrelated 
with the bin radius r,. Here we have applied a variant of 
this test described in detail in BMG97. The field is divided 
into rings of fixed width A6' and distance 9i from the QSO, 
and the variables to be correlated are w{9i), the value of 
the empirical angular correlation in the i— bin, and 8i. The 
results of several binnings are averaged to reduce the depen- 
dence of the significance on the particular choice of AO. Ex- 
tensive Montecarlo simulations show the robustness of this 
approach. 

The main advantage of this test is that it does not rely 
on any particular shape of the correlation function Wqg . It 
just tells us whether the distribution of galaxies is correlated 
with the positions of the QSOs. From it we find that the IJy 
sample shows a positive correlation with galaxies at a 99.0% 
confidence level. Note that the IJy 'North', 6 > sample 
which is also practically complete and does not overlap with 
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log 

Figure 4. Cumulative number counts-flux distribution for the 
HJ and IJy catalogs (we include all the z > 0.3 quasars in the 
original Half-Jansky and IJy catalogs). The solid lines are least 
squares fits 

the sample considered here, was found to be correlated with 
foreground APM galaxies at a 91.4% level (99.1% for the 
red galaxy sample Benftez & Martmez-Gonzalez 1995; see 
also Norman & Williams 1999). A naive but conservative 
combination of both results yields p > 99.9 for the existence 
of positive correlations between the whole IJy catalog and 
foreground galaxies. The existence of correlation for the HJ 
sample is detected at a smaller significance level, p = 84.1%. 

To compare among the different samples wc found more 
convenient the weighted average test of Bartsch et al. (1997). 
This test is based on the estimator rg, 

1 ^ 

^3 = J^^WiaiOi) (6) 

where A'^ is the total number of galaxies and 9i are the 
quasar-galaxy distances. Bartsch et al. (1997) showed that 
rg is optimized to distinguish between a random distribu- 
tion of galaxies around the quasars and a distribution which 
follows an assumed Wqg. 

We will slightly modify the procedure of Bartsch et al. 
(1997) and use a normalized value of rg, which we define as 
r-w = Tg/Inigg where 

= — / w,g{e)ede (7) 

Values of r^, > 1 will indicate a positive correlation, < 1 
negative, and rw = Q absence of correlation. If the galaxies 
are exactly distributed following Wqg, then rg happens to 
be the Monte-Carlo integral of w^g, i.e., rg = 1^2 . Unlike 
Bartsch et al. (1997) we do not merge all the fields into one 
single superfield, but calculate r^ for each quasar and then 
found an average < r^, > over all the sample. This avoids 



giving more weight to a field just because its galaxy density 
is higher, which could bo due to star contamination. 

We set Wqg oc as predicted by Bom'tez & Sanz 

(1999), and proportional to the galaxy angular correlation 
of the APM galaxies (Maddox et al. 1990). This test is very 
robust in the sense that it is only sensitive to changes in the 
shape of Wqg{9); if we multiply this function by am arbitrary 
amplitude, the significance will not change. The values of rw 
together with their significance level Pw , are listed in Table 1. 
The significance is established with 10000 sets of simulated 
fields each with the same number of randomly distributed 
galaxies as the real fields. 

3.1 Is there an 'identification bias'? 

Prom Table 1 we see that the VC sample, the one more simi- 
lar to the PKS sample used in BMG97, has the higher signif- 
icance level, pn, = 99.55% (the PKS itself has pn, = 99.78% 
using this test). However, the union of the VC, HJ and IJy 
samples, which includes only 8 more HJ quasars has a con- 
siderable lower value of rw and of pw If we analyze sepa- 
rately these 8 "non VC" HJ quasars we see that they have 
rw = 0.97132 and pw = 6.72%, i.e. they are strongly an- 
ticorrelated with foreground galaxies at the 93.28% level. 
However, if we look at the other seven HJ faint objects 
which have similar optical and radio fluxes {Bj > 18.6, 
5*2.7 < 0.8) but were included in the VC catalog, they have 
< r^ >~ 1.03805 and p„ = 93.43%, that is, an excess much 
stronger than that detected for the full HJ sample. If we 
put together these 15 faint quasars, we get r^ = 0.99953, 
that is, practically no correlation whatsoever, positive or 
negative. Since magnitudes, radio-fluxes, absolute luminosi- 
ties, redshifts and galactic latitudes are very similar for both 
faint "minisamples" , the only remaining difference seems to 
be their inclusion in the '96 version of the Veron-Cetty & 
Veron catalog, something which only depends on the date 
when the quasars were identified. That happened in 1983 for 
the 7 VC quasars, whereas the "non-VC" redshifts were first 
published in Drinkwater et al. (1997). If we compare the val- 
ues of rw of the 8 "non-VC" objects with the rest of quasars 
(86) from the HJ sample, a Kolmogorov-Smirnov test shows 
that their distributions are incompatible at the 99.31% 
level. This is rather puzzling, and barring a statistical fluke, 
the only remaining possibility seems to be something which 
might be called "identification bias": apparently the first 
radio sources to be spectroscopically identified in a catalog 
tend to be those in regions of higher galaxy density. It is 
possible that due to the positional uncertainty of the radio 
identifications, there was a tendency to start the identifica- 
tion of a radio catalog with those fields which have more 
"candidates" close to the radio position. Although in our 
case the samples are small, the differences among them are 
so significant that one must conclude that any results about 
quasar-galaxy associations obtained with incomplete cata^ 
logs should be considered with great caution. 

The results of BMG97 are a good example: if we look 
at the 86 quasars in the HJ sample considered in that paper 
(where the 8 "non-VC" objects were not included) one ob- 
tains rw = 1.01554 and pw = 99.14%. The total HJ sample 
gives rw = 1.01178 and pw = 97.33%, as we see from Table 
1. Although the galaxy excess does not disappear when we 
use the complete sample, the values of rw and pw drop ap- 



© 0000 RAS, MNRAS 000, 000-000 



6 Bemtez, Sanz & Martinez- Gonzalez 



0.05 




0.01 




Figure 5. The product feCp^ for several omega values, flm = 
0.1 (dashed), 0.4 (dotted) ,1 (solid thin Une). The thick line is 
< bCaS >■ 



preciably. If this happens with a sample which was almost 
(86/94 = 91%) complete, results as those of Tyson (1986) 
, Hintzen et al. (1991), or Burbidge et al. (1990), which 
use QSO samples extracted from incomplete compilations 
as those of Hewitt & Burbidge (1980) and the Veron-Cetty 
& Veron (1984) catalogs, may be seriously biased, explain- 
ing the extremely strong amplitudes of Wqg found by these 
authors. 

Therefore, it seems clear that valid samples for quasar- 
galaxy correlation measurements should be complete or ran- 
domly extracted (e.g., depending on the celestial coordi- 
nates) from complete catalogs, as the IJy and HJ. Note that 
the IJy has the highest value of for the samples in Ta- 
ble 1, i.e. it is the most strongly correlated, in qualitative 
agreement with the magnification bias effect, although pn, is 
not very high due to the smaller size of the sample, which 
increases the statistical uncertainty. 



4 COMPARISON WITH THEORY 

To compare with theory, we shall assume that the biasing 
factor b = Sg/Sm is approximately constant within the rele- 
vant angular range. Two estimates of Wqg will be considered, 
that of Sanz, Martinez-Gonzalez & Bem'tez (1997), which 
takes into account the nonlinear evolution of the power spec- 
trum using the ansatz of Peacock & Dodds (1996), and the 
more recent estimation of Bemtez & Sanz (1999). 



4.1 Power spectrum-based estimate 

In Sanz, Martmez-Gonzalez & Bem'tez (1997) we showed 
that the the value of C^f in Eq. (^ practically does not 



Figure 6. Quasar-galaxy correlation for the HJ (to the left) and 
IJy (to the right) samples. The error bars are poissonian and the 
bin width is approximatel y 1 '. The dotted line corresponds to 
(oe — 1) < bCfxS > (see sec j4.l| ). The two dashed lines correspond 
to the estimation of Wqg of Sec. [4.2| with Q/b = 0.3 (lower) and 
f2/6 = 1 (upper). We can see that in both cases the observed 
correlation is much larger than the expected one. 



change if, instead of using the real galaxy and quasar dis- 
tributions for the calculations, it was assumed that all the 
galaxies and quasars were respectively placed at redshifts zi, 
and z f . Our calculations are normalized to the cluster abun- 
dance (Viana & Liddle 1996) as = o.en-^ ^^+O-^s^-o isn^^ 

For consistency we have to normalize the large-scale bi- 
asing factor as 6 = fg"^- '^^^ above, ujgg oc bC^s, so this 
product contains all the information relevant to the depen- 
dence ofujqg on SI. The variation of b and of the amplitude of 
Cfj,s with n somehow cancel and the product bC^s depends 
weakly on fl, as it is shown in Fig 5, where we plot bC^s 
calculated using delta-functions for the galaxy and quasar 
distributions with Zg = 0.15, Zq = 1.4. We see that there are 
two different regimes for Wqg, one at small scales 9 < 1' and 
the other at larger scales, 9 > 2'. 

The dependence of Wgg on Q, at small scales can be 
approximated as 

Wqg ^ As{n) l0g9 + Bs 

This fit is valid for the angular scales 0.1' < 9 < 1'. As{0.) = 
-0.029r2"'^'''°*5" and Bs (which is equal to the value of Wgg 
at 1') is practically a constant, Bs = 0.029±0.02. As quickly 
drops as Q grows, although the amplitude of the variation is 
relatively small. As = 2.34 at f2 = 0.1 and As = 1 for a flat 
universe. At larger scales, 9 >> 1', Wgg is reasonably well 
approximated by a law of the form 

Wqg Al9~^ 

The amplitude Ai also weakly depends on Cl, with Ai = 
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0.0290}-''^, but unlike the situation on small scales, it grows 
with the value of SI. 

Since the variation of Wqg with SI is small, the obser- 
vations will only be compared with the product (oe — 1) < 
6C,i5 >, where < bC^s >= dflbCfxs{^)- This comparison 
is shown in Fig 6a, where we plot Wqg as measured from the 
HJ and IJy samples. It is clear that the observed correlation 
is much stronger than the model prediction. 

To estimate this mismatch we perform a maximum like- 
lihood fit to the data of the form 

wgg = A{ae - 1) < bC^6 > (8) 

and leave A as a, free parameter. We found that for the 
HJ sample A = 12.7 ± 6.5, where the error limits enclose 
a 68% confidence interval, and A > 1 at the 98.40% level. 
The IJy sample also displays a much stronger correlation 
than expected: A = 10.2 ± 4.5, and A > 1 at the 99.08% 
level. These results are not appreciably changed by varying 
bCi_is{fl) within the range 0.1 < Q < 1. 



4.2 Model independent estimate 

Benftez & Sanz (1999) propose a model- free estimation of 
Wqg, summarized in Eq. (j^. We can only apply it approxi- 
mately, since it assumes that the foreground and background 
galaxies are localized in 'thin' slices, with the angular cosmo- 
logical distance practically constant. This is not strictly the 
case here, but a rough estimate can be found as Wqg oc Wgg 
and Aqg « AggQ{zf, zi,)Azf{a — 1)^, where for Zf, and z/ 
we will take the median values of the QSO (zj = 1.4) and 
galaxy populations (z/ — 0.15) respectively. Also for the 
galaxies Az/ ~ 0.2, Agg ~ 0.44. Since the proportionality 
factor is Q(0.15, 1.4) ^ 0.45 (Bem'tez & Sanz 1999), and we 

A o 

expect to measure a correlation amplitude of ~ 0.04^. 
A maximum likelihood fit to the data using a function of 
the form Wqg oc A9 ' , taking a fixed 7 = 0.688 and leav- 
ing A as a free parameter, yields ^^"ti = 0.39 ± 0.20 and 

= 0.31 ± 0.14, again much higher than expected. 
Given the redshifts of the background and foreground pop- 
ulations, including a cosmological constant A would not 
change significantly the results. 



5 DISCUSSION 

Although the observed results seem to qualitatively agree 
with the lensing hypothesis (for instance, the amplitude of 
Wqg divided by ae — 1 is almost constant for both the IJy 
and the HJ catalogs), they cannot be quantitatively recon- 
ciled with it. For instance, from the results of the above 
paragraph and the estimate of Bemtez & Sanz (1999) one 
would need ^ ~ 10. There seems to be plenty of evidence 
for f2 ~ 1/3 (see e.g. Bahcall et al. 1999), and therefore the 
above ratio would imply absurdly low values of the biasing 
parameter b. In spite of some hints of 'antibiasing' in low f2 
LSS simulations (Jenkins et al. 1998), most measurements 
point at values of fe ~ 1 for the APM galaxies (e.g. Frie- 
man & Gaztafiaga 1999), although on scales several times 
larger than the ones considered here (i? < 1.5h~^ Mpc at 
z < 0.15). 



It seems difficult to believe that the weak lensing the- 
ory is wrong by an order of magnitude. For instance, the 
calculation by Williams (1999) taking into account second- 
order effects only changes in ~ 10% the first-order results. 
Are there other possible explanations for these high values 
of the QSO-galaxy correlation? 

One option often considered is absorption by Galactic 
dust (Norman & Williams 1999). It is easy to see that this 
effect would lead to a positive correlation between galaxies 
and high-z quasars: Let's suppose that 5t{6) represents a 
(small) fluctuation of the galactic extinction in the sky. The 
induced fluctuations in the galaxy and quasar number den- 
sity, Sug and Snq respectively would be Suq ~ —aqSr and 
Sug « —agSr, where a, and ag are the slopes of the opti- 
cal number counts magnitude distributions of quasars and 
galaxies. Therefore, 

Wqg =< SnqSUg >« aqagCrr (9) 

where Ctt is the dust-dust correlation function in e.g. 
Schlegel, Finkbeiner, & Davis (1998). Although the quasar- 
galaxy correlation Wqg also increases with the value of the 
slope aq, there are significant differences with the correla- 
tion generated by the magnification bias effect: Wqg is always 
positive, even for values of a, smaller than 1. This could al- 
low to distinguish between the contribution of each effect to 
the total correlation. The best observational results on Ctt 
available so far, those of Finkbeiner, Schlegel & Davis (1998) 
do not resolve adequately the angular scales on which Wqg 
is measured in this paper. A tentative extrapolation of the 
Ctt presented in the above mentioned paper is too fiat to 
explain the measured correlations. Future measurements of 
Ctt on small scales will serve to establish the contribution 
of TOgg to the total, observed Wgg. 

It is also interesting to explore the contribution to Wqg 
on relatively large scales, > 1', from strong lensing effects. 
Let's suppose that a fraction fs of the quasars in a sam- 
ple is strongly lensed (not necessarily multiply imaged) by 
individual foreground galaxies. These galaxies, which will 
be very close to the quasar positions, are correlated with 
other galaxies, following a correlation function Wgg. If the 
typical quasar-lens distance is considerably smaller that the 
galaxy-galaxy correlation scale, this would indirectly cause 
a "strong- lensing" quasar-galaxy correlation Wqg ~ fsWgg. 
Since the "weak lensing" correlation estimated above is 
uij^ « 0.02 — 0.04u;gg, relatively small values of fs could 
generate a correlation Wqg of comparable amplitude. A quick 
inspection of the Digital Sky Survey images around the 
quasars in the HJ sample reveals several cases of very close, 
6 < few arcsec 'associations' in a sample of ~ 100 quasars. 
Unfortunately, and due to the poor quality of these images 
(some of the 'associated' objects may be stars or defects), it 
is difficult to give an exact value of fs for this sample. How- 
ever, if an appreciable fraction of these objects are galaxies, 
this effect would contribute significantly to the total am- 
plitude of Wqg. Obviously this point merits a more detailed 
consideration with improved observations. 

Last but not least, another possibility is that the ob- 
served correlation Wqg is affected by unforeseen systematic 
effects, or due to a statistical fluke. This sounds very un- 
likely given the variety of quasar and galaxy catalogs con- 
sidered in Wqg studies (Bem'tez 1997). But it should not be 
forgotten that many of the positive results reported in the 
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literature are based on bright, radio-loud quasars, and they 
are so scarce in the sky that overlaps among the samples are 
unavoidable. 

As in most scientific controversies, the solution will 
eventually come from more and better data: ongoing sur- 
veys as the Sloan Digital Sky Survey (Gunn & Weinberg 
1995) or the 2dF (Folkes et al. 1999) will provide huge QSO 
and galaxy catalogs, from which it will be possible to select 
complete samples well suited to measure Wqg. 



6 CONCLUSIONS 

Gravitational lensing predicts an enhancement of the den- 
sity of background QSOs around foreground galaxies. We 
measure this QSO — galaxy correlation Wqg for two complete 
samples of radio-loud quasars, the southern IJy and Half- 
Jansky samples. The existence of a positive correlation be- 
tween distant z ~ 1 quasars and z ~ 0.15 galaxies is con- 
firmed at a p = 99.0% significance level (> 99.9% if previous 
measurements on the northern hemisphere are included). A 
comparison with the results obtained for incomplete quasar 
catalogs (e.g. the Veron-Cetty and Veron compilation) sug- 
gests the existence of an 'identification bias', which spuri- 
ously increases the estimated amplitude of Wqg for incom- 
plete samples. This effect could explain most of the very 
strong quasar-galaxy statistical associations found in the 
literature. Nevertheless, the value of Wqg that we measure in 
our complete catalogs is still considerably higher than the 
predictions from weak gravitational lensing theory. Includ- 
ing the efi'ects of strong lensing could help to explain this 
discrepancy. 
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